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ABSTRACT 

Using the factorization scheme for the nonleptonic D VVq weak amph- 
tudes, we classify all diagrams which arise m. D ^ V'-f decays and calculate 
them with the help of the hybrid model which combines the heavy quark 
effective theory and the chiral Lagrangian approach. Thus we determine the 
long distance contribution to the amplitudes of Cabibbo allowed and Cabibbo 
suppressed D decays. The calculation of the expected range of the 

branching ratios of nine different D — ^ channels is compared with results 
of other approaches. The present work establishes an increase of the parity 
violating contribution in these decays in comparison with previous analyses. 



1 Introduction 

The study of nonleptonic decays of charm mesons has been a subject 
of high priority for more than two decades and has resuhed in a weahh of 
experimental data, which continues to expand at a remarkable rate. This 
rich source of information plays a decisive role in the development of the 
theoretical treatment of these processes, which are governed by the interplay 
of the weak and strong interactions. 

On the other hand, there is very little information available on the sector 
of flavour changing radiative decays of charm mesons, in which the electro- 
magnetic interactions is also operative. 

Only some preliminary upper limits, at the 10"'* range, for branching ra- 
tios of Cabibbo forbidden decays of type V^'j have been reported so far 
[0]. However, as a result of ongoing efforts p[ it is reasonable to expect that 
the experimental data on D ^ V'j decays will be forthcoming during the 
next few years. The theoretical treatment of these decays must address firstly 
the question of the relative importance of short and long distance contribu- 
tions. For the decays studied here, the short distance process of relevance, 
the c — i> M7 transition which is driven by the magnetic penguin diagram, 
is exceedingly small being suppressed by GIM cancellation and small CKM 
matrix elements [Q. The inclusion of gluonic corrections 0, |^ increases the 
free quark c —>■ wy amplitude by several orders of magnitude. However, even 
after taking this increase into consideration the inclusive branching ratio due 
to the c ^ wy short distance penguin reaches only the 10~^ region, which 
is still much smaller than the effect of long distance contributions. Thus, in 
order to estimate these decays, one must concentrate on the treatment of the 
long distance dynamics involved in the D —>■ V'j transitions. 

During the last few years several papers have appreared in which the 
D —>■ V'J transitions were considered. In ref. the first comprehensive 
phenomenological study of the various D V-y has been presented, using 
mainly the techniques of pole diagrams and vector meson dominance. Other 
approaches include the use of the quark model picture and of effective La- 
grangian p, , the use of QCD sum rules and the hybrid model approach, 
which combines heavy quark effective theory and chiral Lagrangian p, |^, p!0|] . 
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As it will be emphasized in the last section, the existing predictions of the 
various theoretical attempts are quite divergent for some of the D —>■ V'y 
modes which underlines the urgent need for experimental data as well as for 
the development of a reliable model. 

It should also be mentioned that these decays offer certain opportuni- 



ties to search for signals of physics beyond the standard model P, H, 11 



although some of the proposed tests could be affected by the long distance 



contributions embodied in the c — > M7 transition |]I0 . 

In the present paper we aim for a more systematic and comprehensive 
treatment of these decays than previously undertaken, employing the for- 



malism of the hybrid model [0, |T3[ for treating the D V'y transition 
[^, [1^. In Section 2, we present the details of our approach and we define 
the approximation used, in Section 3 we give the explicit form of the am- 
plitudes and we conclude in Section 4 with a discussion and a comparative 
presentation of our numerical predictions. 

2 Model description 

We treat the radiative decays D V'y as originating from the non- 
leptonic transition D VVq, followed by the conversion Vq — > 7 via the 
vector meson dominance mechanism. Although a similar scheme has been 
considered also in previous papers 0, H, |1^ there is no systematic treatment 
which includes all possible diagrams within the chosen approximation. In 
the present paper we adopt the factorization approach for the D VVq am- 



plitude, following the formalism advanced in for the nonleptonic decays 
of D, Dg mesons (BSW scheme). We shall not repeat here the arguments for 
using the factorization approximation, as these have been amply discussed 



in the literature (see, e.g. [|1J, |T5[). We are aware that the nonfactorizable 



contributions could also play a role as it may be the case for certain D |]T6 



and B nonleptonic decays. However, at the present time, before any ac- 



tual measurements of D — > V'y exist, we prefer to limit ourselves to a simple 
scheme, and to keep our approach as transparent as possible, awaiting the 
confrontation with experiment. The factorization amplitude for D VVq in 



2 



BSW scheme is calculated using the effective nonleptonic weak Lagrangian 



Gf 

= --^Vugy*g^ [ai{uqi)''{qjc)f, + a2iuc)^{qjqiy'], (1) 

where {ipiip2)'^ = ~ 7^)V'2, qi,j represent the fields of or s quarks, 

Vij are the CKM matrix elements and Gp is the Fermi constant. In our 
calculation we use ai = 1.26 and 02 = —0.55 as found in [|14| . 



In Eq. (|I]) the quark bilinears are treated as interpolating fields for the 
appropriate mesons. In order to calculate the matrix elements we use as 
before p], [1^] the hybrid model which combines the heavy quark effective 



and chiral perturbation theory ||I2|, |T^. The relevant hadronic degrees of 



freedom within this framework are the charm pseudoscalar {D) and vector 
(D*) mesons and the light pseudoscalar (P) and vector (V) mesons. In the 
factorization scheme (" vacuum insertion") which we use, the D Wo 
amplitude is schematically approximated as follows 



{VVomq,nqkc),\D) = {V\{q,q,r\0){Vo\{qkc),\D) 

+ {Vomq,r\0){V\{q,c),\D) 

+ {vvomqjr\Q){Qmc),\D) , (2) 

where the first two terms are the spectator contributions, in the following 
denoted by Aspecy and Aspec,v, respectively, and the third term is the weak 
annihilation contribution, denoted by AAnnih- 

In the three terms of Eq. (||), the Vq meson oj and $) produced in 
the transition is allowed to convert into a photon through the vector me- 
son dominance (VDM). The diagrams thus contributing to the amplitudes 
Aspec,'yi Aspecy ciud AAnnih are shown in the Figs, la, lb and Ic, respec- 
tively. However, as a result of the specific form of the strong Lagrangian of 
the heavy particles (see Eqs. (11), (19) of JlOl)) there is also direct emission 
of the photon from the initial D meson, as exhibited in diagrams (C) and 
(D) of Fig. la. 

The square in each diagram of Fig. 1 denotes the weak transition due 
to the effective Lagrangian (|I|). This Lagrangian contains a product of 
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two left handed quark currents {ipiip2)^i each denoted by a dot on Fig. 1. 
In our model, the left handed currents will be expressed in terms of the 
relevant hadronic degrees of freedom: D*, P and V. In our notation 
in the diagram [B), for example, the hadronic current J2 creates V meson, 
while the hadronic current Ji annihilates D and creates Vq at the same time. 

In ref. |T^ (FS) which also uses the hybrid model for these decays, most 
diagrams exhibited in Fig. 1 have already been calculated. We shall not 
repeat this calculation here and shall combine the results of FS with those 
of our systematic approach, to obtain the full D V'y amplitude in the 
factorization approximation. We rely on FS as a complementary source for 
various basic expressions giving here only those formulae which are directly 
necessary for the calculations of the present paper. 

The principal contribution missing in FS is due to diagram (B) of Fig. 
la. As it turns out, the inclusion of this parity - violating (PV) contribution, 
alters considerably the numerical values of the FS amplitudes in the PV 
sector and leads to the set of predictions for these decays exhibited in Table 
2, which will be discussed in the last Section. The relevant expressions needed 
for diagram (B) are given below. 

The weak current Ji of diagram (B) (see Fig. 1) annihilates quark c and 
creates one of the light quarks q {u, d or s): Jf = q^^{l —7^)0. Under chiral 
SU{?>)l X SU{?>)r this quark current transform as (3^, At the hadronic 
level we impose the same chiral transformation and we require the current 



to be linear in the heavy meson fields D"^ and D*^ Pq , |T2|, |13 



(3) 



where a = fH^/rfiH and ai and 02 are free parameters, which have to be 
determined from the experiment. The current (|^) is the most general one in 
the leading l/m^ order of HQET and next to leading order of chiral pertur- 
bation theory [jl8|, |12|, |13]. Here both the heavy pseudoscalar and the heavy 
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vector mesons were incorporated in a 4 x 4 matrix Ha 



i+MPLr-Pal,) 



(4) 



where a = 1, 2, 3 is the SU{'i)v index of the hght flavours, and P*^, P^, 
annihilate a spin 1 and spin heavy meson Qqa of velocity f , respectively. 
The fields V and u incorporate the light pseudoscalars and are given in FS. 
The field p incorporates the light vector mesons 



■ 9v 



i 



P 



y2 



(5) 



where gy = 5.9 was fixed in the case of exact flavour symmetry (see e.g. 
references in FS). 

The weak current J2 of diagram (P) (see Fig. 1) creates the final V 
meson. Its matrix elements are given by [H 



(6) 



{V{ey,q)\.h,\0) = e:iq)gyiq' 



where the couplings gy^rriy 
mesons: gpijn^p) 



are measured in the leptonic decays of the 
(7p(0) = 0.17 GeV^ gUml) ~ g^O) = 0.15 GeV^, 
g^{ml) ~ g^{0) = 0.24 GeV^ and gx' = {,rnK'/mp)gp. 

For the calculation of the amplitude, we also need the 7 — Vq interaction 
Lagrangian, given by the vector meson dominance like in P, pX[ 



e 

72 



Cy^ -- 

where P^ is the photon field 



(7) 



3 Decay amplitudes 

In order to facilitate the incorporation of the results of FS we shall adopt 
here their notation of the amplitudes, namely A^pQ, A^py for the parity - 
conserving and parity - violating parts, where i denotes classes of diagrams 
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identified below. The general gauge invariant amplitude D{p) — > V{pv)'y{q) 



IS 



G 

A{D{p) ^V{e^v),PiV)h{eiy),q)) = e^Vuqy*q^{e^,^o.pq^e*(^^p''e\^^Apc 

+W{V) ■ -Viv)) - ij>(v) ■ g)(e(y) " ^\^))]Apv}, (8) 

with Apc = A^pc + A^ji(j + A% and Apv = A^py + A^jiy + y4|/^. Now, con- 
cerning the classification of diagrams, Apc* ^"^^ denote the contribution from 
diagrams {A) and (C) (see Fig. 1) which encompass the D* — > transition, 
while A^p(j will denote diagram (G) which contains the P V'j transition. 
Ap^, Apy denote the contribution of the long distance penguins described 
by Fig. 1 of FS. On the hadronic level they are represented by the diagrams 
{E) and (F), respectively, and they vanish in the exact SU(3) fiavour limit as 
shown in FS. We also define two other classes of parity violating diagrams: 
Apy will include the bremsstrahlung-like diagrams (D) and (H), whereby 
the photon emission is due to the direct coupling to charged initial D or final 
V mesons. Finally, Apy, which has not been studied before, will denote the 
contribution represented by the diagram [B). As this contribution will be 
considered in detail below, we present its quark level picture in Fig. 2: Fig. 
2a represents the diagram (B) for the decays proportional to ai, for exam- 
ple while Fig. 2b represents the decays proportional to 02, for 
example Dq —>■ K*^'~f. 

Now we turn to the calculation of the amplitude for the diagram (B), 
while the remaining contributions have been studied in FS. First we parametrize 
the matrix element (Vo| Jfl-D) in the usual way [|TU| 



Vo{q,evoM\D{p))= (9) 
mp) + myQ 



+ z(mD + myo) 



where Q = p — q. In fact the diagram {B) contributes only to form factors Aq, 
Ai, A2 and ^3, which will be determined using the model described above. 
The form factor V{Q'^) gets contribution from the diagrams [A) and (C) of 
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Fig.l, its results are explained in FS, therefore we will leave this contribution 
aside here. In order that these matrix elements be finite at = 0, the form 



factors satisfy the relation |19 



MQ') - ^^^^-^MQ^) + "^^^^f^Mm = , (10) 



and yl3(0) = Ao(0). 



Using the currents and and the relation (10) we determine the 



amplitude of the diagram {B) for ^ p+$ as an example: 

G 

^pv iD+{p) p+(p(y), e(y))$(g, €($))) = -^aiV;^Vud{mD + m$) 
X Aim^y) - J^^^^iP + ^rMm'v)) 9v e\y,,. (11) 

According to the vector meson dominance (|^), the amplitude for 
is obtained, if the polarization e*^^ is repaced by e*^^^eg^ / {?>m%) . However, the 
amplitude for p+7 decay should satisfy the gauge invariance condition. 

It was found |2ll for the case of i? — > K*'y decay, that it is useful to 
analyze the heavy meson decays into VVq in terms of helicity amplitudes 

of the two final vector meson: A and ^00 • Thus, the application 

of gauge invariance condition to D VVq decay, with Vq — 7 conversion, 
means that the Aoq{D VVq) helicity amplitude must be discarded. Under 
a gauge transformation as implemented by ej^^-, g^, we derive the following 
general condition for the D — > VVq — >• V'j decays 



mjj — m? 



[mr, + my^y 



(12) 



imposed for the decays of type (_B) graphs, presented on Fig. la. Conse- 
quently, the amplitude can be expressed in terms of the form factor 
Ax{m\) only. 

Now we determine the form factors y4i(my) for the diagram [S) using 
the current (0) and parametrizing it in the form of (^) . The weak current (|^) 
determines the form factor in the heavy quark limit, i.e. at the maximum 
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momentum transfer Q 



2 

max 



{tud - mvof 



max I 



V2 mo + myo 



(13) 



We assume the pole dominance behaviour of the form factors [[12], [T^ and at 



rriy we determine 



-gvV2ai 



D 



(mu~mvo? 
°i+ 



+ myo 1 - 



(my)^ 
n 



(14) 



where Di+ is the mass of the qc = 1+ bound state. We use the masses 
of sc and dc bound states to be 2.53 GeV and 2.42 GeV as in [|13[. The 
free parameter ai is determined by using the average of experimental Ai(0) 
values for D+ Mvi and D+ K*Hvi. We obtain = 0.171 GeV^''^ 
and use this value for the prediction of all D V'y decay rates. 

Using the formalism described above, with Eqs. ([TTl) , ([121), (|l^) , we obtain 
Apy for the various decay channels. 

The Cabibbo allowed decay amplitudes, which are proportional to the 
product |V^dV7g|, are: 



+ 



mi 



{mo + mp)\A^''{ 



3m,^, 



(m^ + mJ|Af-(m^.)l 



m 



D 



m 



K* 



(15) 



ai 



3md 



■(mz, + m^)|Af=*(m^)| 



■ (16) 



The Cabibbo suppressed amplitudes Apy proportional to the Cabibbo 



factor IKmKsI ^^^6 



-ai 



9l 

m\ 



{mo + mp)\A^\n?^\ 



^(mz, + mJ|Af-(mJ)| 



m 



D 



(17) 
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6ms, 



mjj — m 



K* 



a^Md' - A) 



02 

71 



9l 
rrii 



[mD + mp)\A^''{ml)\ 



+ 



3m2 



m|) — w? 



(19) 



02 

V2 



■{mD + rnp)\A^^{ml)\ 



(20) 



>1^pVP^^7) = -«2 



{mo + ■mp)\A^''{'m. 



+ 



(mD + ma;)|>lf'^(m| 



(21) 



For completeness, we give also the parity violating parts of the amplitudes 
for doubly suppressed decays proportional to iV^usV^^I- 



X*+7) = 



3ml 



'-^{mn + mp)\A?^{ml. 
mi 



■{mo + m^)\A^^' {mj^.)\ 



mjj — m 



-,(22) 



K* 



i^*°7) = 02 



9p9K 



mt 



-{mD + mp)\ A^^ {m\, ) \ 



+ 



9lo9k' 
3ml 



{mD + m^)\A'^'^{mii, 



m% — m\. 



■ (23) 
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4 Discussion 



We present numerical results for the amplitudes Ap^, Ap^, -A^pci Apy, 
PY and jA.^pY 



Apv and Apv in Table 01, where A^ denotes 



■A'pciv) = e^K,,K;.^PC(y) ■ (24) 



and i runs over the nine decays studied. The amplitudes ApQ are calculated 
from Eqs. (15) - (23), while for the rest of the amplitudes we use the values 
obtained in FS. Although we have the values of Table 1, which encompass all 
the amplitudes arising from our factorization scheme and hybrid model, we 
cannot predict at this stage definite values for the decay rates. This is due to 
the fact that the signs of several constants entering the expressions of Ap^, 
Apy (like A', A, Cyvu, gv, ca, defined in FS and here) are not determined yet 
from the experimental data. Thus, we must contain ourselves in the present 
only to the range of values, which are obtained by assuming all possible 
relative signs for the various constants. We also remark that Ap^, Apy are 
usually one to two orders of magnitude smaller than the other amplitudes 
and will affect the rates only in the case of cancellations occuring among the 
rest of amplitudes. 

We calculate the branching ratios of the D V'y decays with the help of 



r(D ^vj) = ^ 



2m 



D 



{\Apc\' + \Apy\'), (25) 



and we present the possible range of values for the branching ratios in Table 
2. We compare the present results, denoted by (a), with the results obtained 



in previous approaches. The results of the approach |]T0[ are denoted by (b), 
the results of |Q by (c), and the results of by (d) (we specify Oi = 1.26 and 
02 = —0.55 in their formulas). The quark model calculation of predicts the 
branching ratio for the Cabibbo allowed decays BR{D^ K*^'-/) = 8.6- 10~^ 
and BR{Df — > p'^'y) = 2.1 ■ 10~^, which are smaller in comparison with the 
results we obtain here, as well as compared with the predictions of H] and 
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1^. On the other hand, the calculation of 0, which also uses a quark model, 
leads to a larger branching ratio BR{D^ —>■ K*^'-/) = 1.1 ■ 10~^, which is an 
order of magnitude larger than obtained in ||^ and closer to our estimate. 
We notice that the parity violating amphtudes calculated within the present 
approach have changed significantly in comparison with the results of FS. 
Overall, inspection of Table 2, indicates that our predictions in the present 
paper are closest to those of When measurements of a few channels will 
be available it will be possible to adjust the range of the various branching 
ratios and to make firmer predictions, hopefully allowing to select the best 
suited model. 

The measurable ratio Apc/Apv can also be used to distinguish between 
the various models. However, since at the present stage we cannot specify 
the relative signs of the various components of each of these amplitudes, we 
are unable to make any sensible statement about these ratios. 

We summarize our results as follows: we have presented a calculation of 
the radiative D —>■ Vj decays using a model which contains all classes of 
diagrams arising from the factorization approach for the D VV'^ ampli- 
tude, from which the radiative decays are obtained by use of vector meson 
dominance. In the calculations of the various matrix elements, we use a hy- 
brid model [12, 13] which combines heavy quark techniques with the chiral 
Lagrangian. In view of uncertainties related to the coupling constants in- 
volved, we can predict at this stage only ranges for the branching ratios of 
the various decay channels, which are given in Table 2. We emphasize that 
the Cabibbo allowed decays D° — > K*^'j and — >• p+7 are calculated to be 
fairly frequent, with branching ratios of a few times 10~^ and we expect their 
detection soon. Some of the Cabibbo suppressed modes, like 1^+'° p+'°7 
may also occur with branching ratios close to 10~^. Experimental results on 
these modes are eagerly awaited and will certainly contribute to clarify the 
long distance dynamics leading to these radiative decays. 
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Figure Captions 

Fig. 1. Skeleton diagrams of various contributions to the long distance 
decay D V^y. The spectator diagrams of type Aspec,-y (see Eq. (Q) in 
the text) are shown in Fig. la, the spectator diagrams of type Aspec,v are 
shown in Fig. lb and the weak annihilation diagrams AAnnih are shown in 
Fig. Ic. The square in each diagram denotes the weak transition due to the 
effective Lagrangian Cy^ (|I|). This Lagrangian contains a product of two left 
handed quark currents {'ipiip2)'^, each denoted by a dot. Different diagrams 
are denoted by (A) - (H). Their contributions to the amplitudes ApQ and 
Apy are specified in the text. 

Fig. 2. The quark level picture of the diagram [B) of Fig. 1 (amplitude 
Apy): Fig. 2a represents the decays proportional to ai, for example — > 
p^7, while Fig. 2b represents the decays proportional to a2, for example 
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^ p+7 


0.59 
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0.02 


1.6 


1.3 


0.03 


^ A^*+7 


0.41 


2.3 


0.02 


2.1 


1.2 


0.04 


^ ir*+7 


0.16 


0.42 





0.43 


0.37 





/}0 ^ j^*0^ 


0.33 


0.32 








0.28 






Table 1: The parity conserving and parity violating amplitudes for charm 
meson decays in units 10^^ GeV^^ . The amplitudes ^pj-^py get contributions 
from different diagrams in Fig. 1: (A) and (C) contribute to J^pc^ (D) and 
(H) contribute to A^py, (G), (B), (E) and (F) contribute to Apc, A%, A^JS 
and Apy, respectively. The amplitudes ApQ, ApQ, ApQ, Apy and Apy 



were calculated in [0], while the amplitude Apy represents the additional 
contribution calculated here. The first two decays are Cabibbo allowed, while 
the last two are doubly Cabibbo suppressed. 
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L> — > V "J 
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(0.1 - 1) 
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DO ^ CJ7 


(0.1-0.9) 


(0.02-0.8) 


~ 0.2 




D" ^ $7 


(0.4-1.9) 


(0.04-1.6) 


(0.1-3.4) 




D+ ^ p+7 


(0.4-6.3) 


(1.8-4.1) 


(2-6) 


0.43 


D+ ^ fs:*+7 


(1.2-5.1) 


(2.1-3.2) 


(0.8-3) 




D+ ^ fs:*+7 


(0.03 - 0.44) 


(0.12-0.25) 


0.1-0.3 






(0.03-0.2) 


(10-5-0.08) 


~ 0.01 





Table 2: The branching ratios for D V'y decays. The first column (a) 
contains the results of the present approach. The next three columns present 
the results of Ref. ^ (b), Ref. g (c), and Ref. (d). The first two decays 
are Cabibbo allowed, while the last two are doubly Cabibbo suppressed. 
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